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Abstract 


A weakly nonlinear theory is developed to study the interaction of 
Tollmien-S chlichtin g (TS) waves and Dean vortices in curved chann el flo w. 
The predictions obtained from the theory agree well with results obtained 
from direct numerical sim ulations of curved ch ann el flow, especially for low- 
amplitude disturban ces. Some discrepancies in the results of a previous theory 
with direct numerical simulations are resolved. 


1. Introduction 

Interactions between Tollmien-Sclilicliting (TS) waves and streamwise vortices can be an impor- 
tant part of the transition process (ref. 1). Kim and Moser (ref. 2) and Singer, Reed, and Ferziger 
(ref. 3) showed that weak streamwise vortices can alter the mode by which the flow undergoes tran- 
sition. Nayfeh and Al-Maaitah (ref. 4) used Floquet. theory to study the influence of Gortler vortices 
on the growth of oblique TS waves. They found that there are two components to the solution; one 
is stabilized and the other is destabilized by the presence of the vortices. Malik and Hussaini (ref. 5) 
studied the interactions of Gortler vortices and TS waves via direct numerical simulations. They 
found qualitative agreement with the results of Nayfeh and Al-Maaitah (ref. 4). 

A series of studies on the interactions of TS waves and streamwise vortices in curved channel 
flows at asymptotically high Reynolds number have been made (refs. 6-9). Daudpota, Hall, and 
Zang (ref. 10) (hereinafter referred to as DHZ) developed a weakly nonlinear interaction theory to 
study the interaction at a finite Reynolds number. They employed a multiple-scale version of the 
approach developed by Stuart (ref. 11) and Watson (ref. 12) to derive two coupled Landau equations 
for the perturbation amplitudes of the streamwise vortices (called Dean vortices in curved channel 
flow) and TS waves. A comparison of their theoretical predictions with the results of direct numerical 
simulation (ref. 13) suggests that their results are in error with respect to the influence of the TS 
wave on the Dean vortex. The resolution of this discrepancy is a major motivation for undertaking 
this work. 

In section 2 we present a slightly different formulation (see refs. 14 and 15 by Herbert) of a weakly 
nonlinear interaction theory and describe how the resulting equations are solved. We also include 
suggestions to guide the intelligent use of the theory. In section 3 we compare the results from the 
current theory with results from DHZ. In addition, direct numerical simulation is used to verify that 
the current theory correctly predicts the overall behavior of the flow. Finally, in section 4 we draw 
conclusions. 

The authors would like to thank Thorwald Herbert of Ohio State University, Q. Isa Daudpota, 
formerly NRC at the Langley Research Center, and Philip Hall of the LIniversity of Manchester,. UK, 
for many useful discussions. 


2. Mathematical Formulation 

2.1. The Basic Equations 

The incompressible Navier-Stokes equations in cylindrical coordinates (r*,0,z*) are written as 
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and p* and v* represent the constant density and kinematic viscosity, respectively. The asterisks 
indicate dimensional quantities. The geometry of the problem of azimuthal flow between infinite 
concentric walls of outer radius r* and inner radius r* is illustrated in figure 1. The wall boundary 
conditions require that 

Ur = u; = U* = 0 (r* = r-,r*) (7) 



A solution to the equations gives (U* ,U * , P*) as 
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Here, U* is the bulk velocity and 
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The quantity r* = (r* + r*)/2 is the centerline radius, h* — (r* — r*)/2 is the channel half-width, 
and the ratio r/ — r*/r * describes the channel curvature. 

We nondimensionalize all spatial coordinates with the channel half-width, velocities with the bulk 
velocity, and pressure with p*U* . The temporal scale is h* /U* . The Reynolds number is defined as 

R=UE 

The nondimensional equations are easily obtained from the dimensional ones by replacing all starred 
quantities with their corresponding nondimensional unstarred ones and noting that p — 1, h — 1, and 
v — 1 / R. The expression 1 — ?/ appears often below and is denoted by A. When A = A c = 2.179xl0 - ' 5 , 
the minimum Reynolds number for instability of TS waves is the same as that for Dean vortices 
(ref. 16). 

We follow closely the perturbation method introduced by Herbert (refs. 14 and 15), extending it 
to the case of interacting disturbances. The steady solution is perturbed such that 
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Substituting the perturbations into the nondimensionalized Navier-Stokes equations, subtracting the 
steady-flow component, and rearranging gives 
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We are considering here the idealized problem in which the solution is periodic in the azimuthal 
(streamwise) and axial (spanwise) directions and the solution evolves in time. In addition, the mean 
pressure gradient is kept constant. 

2.2. Fourier Expansions 

The linearized disturbance equations are obtained by setting the right-hand sides of equa- 
tions ( 14 )— ( 17 ) equal to zero. A solution to the equations for a two-dimensional (2D) TS wave 
can be written as 
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q = qi,ol 

Here K is real and represents the azimuthal wave number, whereas j(t) is a real function defining 
the phase. Any growth or decay of the linearized TS wave is expressed in terms of the amplitude 
such that 

qi,o(r,f) = A(t) qi,o,0,o(f) (21) 
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and ao o is the linear growth rate. Quantities of the form q nm i k are used frequently throughout the 
report. The n and / indices refer to multiples of the TS and/or Dean wave numbers, whereas the m 
and k indices refer to the order of the nonlinear corrections. For instance, qi o 0 0 i s the zeroth-order 
correction of the primary TS wave number, in other words, the TS linear eigenfunction. We define 
the amplitude as 
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where r c is the channel centerline. In general, the amplitude represents the ratio of the magnitude of 
any perturbation quantity and the magnitude of the corresponding linear eigenfunction quantity. The 
amplitude definition can be varied independently of the normalization of the linear eigenfunction. 
This is somewhat more general than that done by Herbert (refs. 14 and 15). The importance of 
this generalization will be made clear later when we discuss alternative definitions of the amplitude. 
Here, we normalize the linear eigenfunction by taking 
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In a similar manner, one can write the solution for the primary Dean vortex mode as 

q = q 0 ,l (M) (25) 


where /3 is the real axial wave number and ((t) — 0 indicates that the linear Dean disturbance is 
stationary. We also write 
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When the nonlinear terms (the right-hand sides of eqs. ( 14)— ( 17) ) are included, the disturbances 
can interact with themselves, with the mean flow, with each other, and with all the relevant 
complex conjugates. This results in the generation of harmonics, mean-flow distortions, and various 
corrections to the fundamental disturbances. Hence, it is natural to expand the perturbation 
variables in the double Fourier series: 
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2.3. Amplitude Expansion 

The nonlinear partial differential equations obtained by substituting .equation (30) into equa- 
tions (14 )— ( 17) are coupled and difficult to solve efficiently. We seek a solution by expanding q n [ 
in the amplitude parameters A(t) and B(t) about the linear solutions given by equations (20) and 
(25). The solutions will reflect the deviations from linear behavior for finite, but sufficiently small, 
amplitudes A and B. In this context, all harmonics of the primary disturbances are considered 
to be forced. Since the linear TS and Dean solutions are 0(A) and 0(B ) , respectively, their first 
harmonics and cross terms will be 0(A 2 ), 0(B 2 ), and O(AB). Higher order harmonics and cross 
terms will be 0(A n B l ); hence, it is reasonable to let 


q rc,/(M) = A^B^ q n j(r,t) (31) 

where all double tilde terms are 0(1) except for those with n — l — 0. Heuristically, one can see that 
the exception with the n — l — 0 term comes about from the fact that the lowest order mean-flow 
distortions are generated from the product of either TS or Dean fundamental disturbances with 
their respective complex conjugates. The fundamental disturbances are O(A) and O(B); hence?, the 
product terms that generate n — l — 0 are 0(A 2 ) and 0(B 2 ), respectively. 

Substitution of the representation in equation (31) into Equation (30) yields 
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We now expand the double tilde representations into sums that are products of ascending powers of 
the amplitude functions with coefficients which are strictly functions of the radial coordinate. In the 
limit, as A — 0 and B — 0, the solutions tend toward the linear results. Only even powers of A and 
B are needed (Drazin and Reid in ref. 17) because of the invariance of the original equations and 
boundary conditions with respect to arbitrary translation in the streamwise and spanwise directions 
and the assumption of periodicity of the solutions in these directions. Hencdj 
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By substitution into equation (32) one obtains the full representation of the perturbation vector: 
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In order that the wall boundary conditions be satisfied at all orders of approximation, we require 
that 


f - n.m.l.k 
Ur 
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at r = r 8 ', r 0 for all combinations of n, to, /, and k. The case n — m — l — k — 0 represents the 
basic flow which has already been considered so that q 0 o 0 0 — 0. 


Consistent expansions of the time derivatives A/ A, B/B, j, and ( are in powers of A 2 and B 2 
so that 
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Equations (36) and (37) are the coupled Landau equations that describe the growth of the 
disturbances. The purpose of the remainder of this section is to determine the Landau coefficients 

2.4. Solution Method 


The solution expansions in equations (32)-(39) are substituted into the perturbation equa- 
tions (14)— (17) . All terms with common exponential factors are grouped together. The simplest 
nonlinear theory is obtained by considering only those terms with Fourier exponents limited by 


— 1 < n and / < 2 


With this restriction, it is only appropriate to consider terms in the amplitude expansions (eq. (33)) 
with to = 0, 1 and k — 0,1. Larger values of to or k give higher order contributions that should 
not be considered at this degree of truncation. In each group with common exponential factors, all 
terms with common powers of A and B are collected. MACSYMA (ref. 18) was used to substantially 
reduce the chance of error in obtaining the final equations. Twelve systems of equations, indexed by 
s, remain to be solved numerically. A description of each of the systems is summarized in table 1. 
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Table 1. Summary of Equation Systems 


Equation system 
number (s) 

Order 

n 

m 

l 

k 

Description 

1 

A 

i 

0 

0 

0 

Linear TS wave 

2 

B 

0 

0 

1 

0 

Linear Dean vortex 

3 

A 2 

0 

1 

0 

0 

Mean-flow distortion from TS 

4 

B 2 

0 

0 

0 

1 

Mean-flow distortion from Dean 

5 

A 2 

2 

0 

0 

0 

TS harmonic 

6 

B 2 

0 

0 

2 

0 

Dean harmonic 

7 

AB 
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0 

1 

0 

TS — Dean cross harmonic 

8 

AB 

1 

0 

-1 

0 

TS — Dean cross harmonic 

9 

A 3 

1 

1 

0 

0 

Self-correction to TS 

10 

B 3 

0 

0 

1 

1 

Self-correction to Dean 

11 

AB 2 

1 

0 

0 

1 

Dean correction to TS 

12 

A 2 B 

0 

1 

1 

0 

TS correction to Dean 


The first two equation systems (s = 1 and 2) are eigenvalue problems of the form 

L n, An, 0,1,0 - [ n («0,0 - *00, o) + * (Vo - *^0,o)] Mq nj0 p 0 (40) 


where 
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The elements of the operator L n ; are given in appendix A. The solution of the eigenvalue problem 
determines the most unstable eigenvalues, aq o — Woo and &o 0 — «7i o 0- The normalizations in 
equations (24) and (29) are used to uniquely define the eigenvectors qi o 0 o( r ) and 4ooio( r )- The 
eigenvalues are the zeroth-order Landau coefficients. 

Equation systems 3 and 4 (i.e., s — 3 and 4) represent the mean-flow distortion caused by the 
TS and Dean disturbances, respectively. For this case, n — l — 0 and the operator L n / is greatly 
simplified. We consider the case where either to or k equals 1, while the other equals 0. The 
mean pressure gradient is constant; hence, p°T°>° = p 0,0,0,1 — q. Analysis indicates that only the 
azimuthal velocity component of the mean-flow distortions is nonzero. Its equation is 
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In similar contexts, Davey and Nguyen (ref. 19) and Herbert (ref. 15) point to an important 
restriction to the theory that results from equation (42). When either of the linear growth rates ao,o 
or &o 0 are negative, the associated homogeneous problem with homogeneous boundary conditions 
can be recast as a Bessel equation through a change of variables. Hence, solutions of the homogeneous 
problem can exist and can contribute an arbitrarily large component to the solution of the forced 
problem. To prevent this difficulty we require both ago and &o 0 t° be nonnegative. This restricts 
the application of the theory to those regimes where both the TS and Dean disturbances are linearly 
unstable. 

The equation systems 5-12 (i.e . , s — 5-12) are all of the form 

[L n j - (2ma 0j o + M a o,o - *«0o,o + 2fc6 0 ,o + |^o,o - *^o,o) ^L]^ln,m,l,k - R s (43) 

where R s is the right-hand-side vector of the s equation system. The right-hand-side vectors consist 
of nonlinear combinations of the solutions obtained at lower order. The elements of the right-hand- 
side vectors R s are given explicitly in appendix B. 

Equation systems 5-8 (i.e., s — 5-8) for the various harmonics require numerically solving a 
sequence of linear, ordinary differential equation systems. In the case of Dean disturbances, we 
found that it is possible that the harmonic of an unstable mode is also in the unstable regime. When 
this happens, it violates the assumption that all the higher harmonics are forced by the primary 
disturbances. Although numerical results can be obtained for such cases, the results are meaningless 
since the harmonic can grow on its own and is not simply forced by the fundamental. We found that 
this circumstance does not occur for the TS waves. 

Equation systems 9, 10, 11, and 12 (i.e., s — 9, 10, 11, and 12) define the corrections to the 
fundamental TS and Dean disturbances and allow us to calculate the Landau coefficients ai o, &o,l ; 
ao i, and b i q, respectively. It is in the solution of these terms that the method developed by Herbert 
(refs. 14 and 15) differs from the standard approach developed by Stuart (ref. 11) and Watson 
(ref. 12). 

The right-hand-side vectors R s for s — 9-12 can be rewritten as 

Rs — R) + K sMq ni o,[,o (44) 


where 


^9 = ai,o — iQ 1,0 
K io — &o,i ~ *^0,1 
K 11 = a o,i - *00,1 
k 12 — &i,o ~ 

The correspondence between the index s and the indices n, to, /, and k is given in table 1. The 
solution vectors q n m \ /ufor these equation systems are the sums of particular solution vectors ; hence, 

q n,m,l,k{r) = K» X° n ,m,l,k( r ) + Xn,m,/,*( r ) ( 45 ) 
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where 


(46) 
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Substitution of equations (44) and (45) into equation (43) leads to 
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Recognizing that 
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when (n — 1, / = 0) or (n — 0, / = 1), one hnds by inspection of equation (47) that 

- (2ma 0 ,o + 2fc6 0 ,o) X° n ,m,l,k = <ln,o,l,o 


or 


0 0 , 1,0 

Xn,m,l,k- 2ma 0)0 + 2kb 0)0 


(48) 

(49) 

(50) 

(51) 


It is straightforward to solve equation (48) numerically to obtain m l k- Substituting equation (51) 
into equation (45) gives 

4„., ,,“(■■) = d, (52) 

In order to determine the Landau coefficients k s , we use the amplitude definitions in equa- 
tions (23) and (28). Together with the expansions in equations (31)— (34) they require that 
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For these to be correct at all orders of approximation, 
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(53) 

(54) 

(55) 
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for all to, k > 0. 

The conditions from equations (55) and (56) can now be used to determine k s . Let the radial 
velocity component of x* m \ ^ be Ur’ m,l ’ k . Equations (55) and (56) are applied in equation (52) to 

. n, 0,1, Of 

(57) 


require 


Hence, 
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The correction to the fundamentals is found by substitution back into equation (45). 


(58) 
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2.5. About the Definition of the Amplitude 

Changing the definitions of the amplitudes in equations (23) and (28) leaves the analysis 
essentially unchanged. Operationally, only the quantity (r c ) /ii™ ’°’ l ’°(r c ) in equation (58) 

needs to be changed. All the numerical computations up to the calculation of k s remain the same as 
long as the normalizations of the linear eigenfunctions in equations (24) and (29) are fixed. Hence, 
one can study the influence of different amplitude definitions on the Landau coefficients without 
repeating the bulk of the computations. (The equations to be solved for higher order Landau 
coefficients depend on the lower order Landau coefficients; so if higher order Landau coefficients are: 
to be calculated, changing the amplitude definition would require recomputing many more solutions.) 
In this subsection, we make some observations regarding how the Landau coefficients vary with the 
amplitude definition and how this information can guide the use of the weakly nonlinear theory. 

All valid definitions of the amplitudes tend toward the same Landau coefficients when the 
linear growth rates ao,o and b o,o tend toward zero. One Can see this by following the approach 
of Watson (ref. 12) and expanding ; k in, equation (48) in powers of 2(?nao o + &&o,o) so that 


L n,m,l,k 


2(ra«o,o + kb 0fi ) 


. 1,-1 


. 1,0 


2(ma 0 ,o 


kbofi )\n,m,l,k 


(59) 


By equating like powers of 2(mao,o + t’i)o,o). we find that at leading order, ; k (r) is proportional 
to the linear eigenfunction q ra oio( r )- Herbert (ref. 15) uses this in his proof that the Landau 
coefficients obtained with his method are the same as those obtained with the Watson method 
(ref. 12) as the neutral curve is approached. Varying the amplitude definition does not change the 

ratio uj'”' m,! ’ i '(r c )/«r’ i:i ’ ! ’ 0 (r c ) in equation (58) if X * m / /,-(?’) is proportional to q raj o,;,o( r )- Hence, for 
2(mao o + kbQ o) — 0, the Landau coefficients become independent of amplitude definition. 

When the quantity 2(raao j o + &&o,o) is not small, x\ m i k( r ) is not approximately proportional 
to q„ 0 ,/,o ( r ) & nd the Landau coefficients can depend strongly on the amplitude definition. Local 
amplitude definitions like the ones used above will generate Landau coefficients that describe the 
local behavior of the disturbance. Locally, parts of the disturbance can grow while other parts 
can decay; hence, the Landau coefficients corresponding to different amplitude definitions may not 
even have the same sign. As noted by Herbert (ref. 15) and others, the dependence of the Landau 
coefficients on the amplitude definition corresponds to a rearranging of the higher order terms in 
the infinite series. Large variations in the Landau coefficients as the amplitude definition is changed 
should warn the user of the theory that the higher order terms are likely to be important for the 
specific problem considered. 

Because integration tends to smooth rapid local changes, amplitude definitions that measure 
global quantities are likely to be more reliable in predicting the overall behavior of the disturbance 
in these cases. Therefore, onf useful alternative amplitude definition is a slight variation of one 
proposed by Herbert (ref. 15). Here, 


A(t) = 


Jr ° q-i,o,o,o(v) • qi,o(v) dr 
/r/q— 1,0 , 0 , 0 ( 7 *.) - qi,0,0,o(?') dr 


(60) 


with a corresponding definition for B(t). Note that q-i^o.ois just the complex conjugate of < 3.1 , 0 , 0 , 0 - 
The expression for the Landau coefficient is then changed from equation (58) to 


K s 


(2ma 0 ,o + 2A,*6 0 ,o) 


Jr° -Xn,m, I,k( r ) dr 

fr ° q - n , 0 ,- l , d r ) • fin, 0 ,L o( 7 *.) dr 


(61) 


Note that the normalization of the linear eigenfunctions need not be changed. 

Because DHZ and much of our early work used local amplitude definitions, comparisons will 
be somewhat simpler if we use a local amplitude definition. Unless otherwise specified, the results 
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reported here used equations (23) and (28) for the amplitudes and equations (24) and (29) for the 
eigenfunction normalizations. 

2.6. Numerical Approach 

All the equation systems are solved using a Chebyshev collocation procedure (ref. 20). The 
left-hand-side operators for all the systems (except the mean-flow corrections) are generated with a 
single subroutine. The operators for the mean-flow corrections need to ensure that only streamwise 
velocity perturbations are nonzero. This was most easily done using a separate subroutine. The 
right-hand-side vectors were sufficiently different from each other that it was more convenient to use 
MACSYMA to write individual subroutines for each of the equation systems. The two eigenproblems 
are solved first with a global solver to obtain the eigenvalue spectrum. The eigenvectors associated 
with the most unstable eigenvalue are obtained with a local iterative procedure. For A < 10 -4 (which 
is the regime studied here), a straightforward evaluation of the mean flow using 64-bit precision leads 
to a substantial roundoff error; hence 128-bit precision was used for the calculation of the mean- 
flow profile. For the stability calculations, we found that 65 collocation points across the channel 
(including the end points) produced Landau coefficients that agreed to at least 3 digits with those 
obtained with 97 collocation points. 


3. Results 

3.1. Predictions of the Weakly Nonlinear Theory 

We present some of our weakly nonlinear results in the context of nonlinear autonomous systems. 
Recall that A measures the amplitude of the TS wave and B that of the Dean disturbance. 
Equations (36) and (37) truncated at third order admit four possible steady-state solutions: the 
trivial solution, finite A with B — 0, finite B with A — 0, and a combined state with finite values of 
both A and B. Each of these states is a critical point and may be classified according to the local 
behavior of the solutions in its vicinity. In the phase plane (B versus A), stable nodes are identified 
as those for which solution trajectories near the point converge to the node, unstable nodes are 
critical points from which solution trajectories diverge, and saddle points are those points for which 
a finite number of trajectories converge to the point while all others diverge. 

Several trajectories are illustrated in figure 2 for a case where R — 6291.67, A = 2.189 X 10 - ' 5 , 
K = 74 257, and / 3 — 4.51. These parameters were used extensively by Singer and Zang (ref. 13) in 
their direct numerical simulations. An explicit fourth-order Runge- Kutta differencing scheme was 
used to trace the subsequent trajectories. The filled circles show the two stable equilibrium points 
at {A — 3.8 X 10 -3 , B — 0) and {A — 0, B — 7.0 X 1 0 — 5 ) , respectively. Trajectories initiated 
near either of these points tend to converge to the respective nodes. The squares represent unstable 
nodes and saddle points. Two initial conditions start close to the saddle point at [A — 2.4 X 10 -4 , 
B — 3.5 X 1 0 — 5 ) and go to different stable equilibrium states. The flow is especially sensitive to small 
changes near the saddle point. The unstable node at the origin is the end point of a semi-infinite 
curve that goes through the saddle point, and separates the region of attraction of the two stable 
nodes. Quite different behavior can be obtained depending upon the initial values of A and B. 

A slightly different scenario exists when we consider R — 5000, A = 2.179 X 1 0 — 5 , K = 10' 5 , 
and /3 — 2.0. In figure 3, pairs of trajectories with initial values of A differing by 0.001 are plotted 
for various initial values of B. In many ways, the qualitative behavior seen here is similar to that 
illustrated in figure 2. The origin is an unstable node, there is a saddle point with nonzero A and 
B , and there is a stable node with nonzero B but A — 0. A semi-infinite curve extends from the 
origin, goes through the saddle point, and separates the domain into a region that is attracted to the 
stable node with [A — 0, B — 7.7 X 1 0 — 4 ) and a region in which A — 0O. This shows how important 
the nonlinear interaction can be. Without a Dean disturbance, the TS wave has unbounded growth; 
however, the inclusion of an additional disturbance at sufficiently large amplitude can completely 
stabilize the otherwise growing TS wave. 
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Figure 3. Trajectories in phase space with R = 5000, A = 2.179 x 10 -5 , N = 10 5 , and /? = 2.0. 

3.2. Comparison With Direct Numerical Simulation 

Here we will compare the results of our weakly nonlinear analysis with those obtained by DNS. 

The numerical simulation code (ref. 13) uses a curved channel variant of the method described 
by Zang and Hussaini (ref. 21) with the nonlinear terms in skew-symmetric form (ref. 22). Table 2 
reports the least stable eigenvalues ( a 0j o — *fifo,o) obtained with a spectral linear stability code and 
the corresponding results of DNS with two different time steps. (Note that in this context K and 
/3 refer to the azimuthal and axial wave numbers of the linear waves.) Sixty-live points across the 
channel are used. Data from the DNS are taken after 100 time steps. Using a time step At — 0.0001, 
the complex growth rates differ from those predicted by linear theory by less than 1 part in 10 4 . In 
the simulations reported below, we are primarily concerned with long-time trends and final steady- 
state solutions; hence, the additional time advancement errors associated with the larger time step, 
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At — 0.01, are not expected to be important. The larger time step gives a maximum Courant (CFL) 
number of 0.01, still well within the numerically stable regime. Additional details of code validation 
studies in both the linear and nonlinear regimes are given by Singer and Zang (ref. 13). 

Table 2. Comparison of Linear Growth Rates 

[All values are multiplied by 10 4 ] 


R 

N 

0 

Linear stability 

DNS with At = 0.01 

DNS with At = 0.0001 

5 000 

10 5 

0 

3.1189 — 4236.91 

3.0736 — 4241.81 

3.1185 - 4236.9i 

5 000 

0 

2.0 

6.6331 + 0* 

6.6336 + 0* 

6.6335 + 0* 

5 000 

10 5 

2.0 

— 138.2126— 16 214.81 

-138.2154- 16 214.81 

-138.2126- 16 214.8* 

10 000 

80 000 

0 

86.4287 — 2709.11 

86.4204 — 2709.41 

86.4291 - 2709.1* 

10 000 

0 

3.0 

21.5591 + 0* 

21.5599 + 0* 

21.5599+ 0* 

10 000 

80 000 

3.0 

-90.6119- 12 993.31 

—90.6131 — 12 993.31 

-90.6119- 12 993.3* 


Sufficient spatial resolution in the simulation can be ensured by using the guideline suggested by 
Krist and Zang (ref. 23). They suggest that “grid refinement is needed in any direction where the tail 
of the energy spectrum reaches 10 -8 of the low-frequency value.” This guarantees that truncation 
errors in the velocity will be less than 0.01 percent. Such detailed resolution is not necessary for 
the purposes of these simulations. We have found that the high wave number end of the energy 
spectrum could be as much as 10 -4 of the low wave number value and still provide results similar to 
those obtained when the Krist and Zang guideline was strictly followed. The relative insensitivity to 
resolution is attributed to the rather minor role that small scales play in the present wave-interaction 
problem. 

We have simulated four cases with R — 6291.67, A = 2.189 X 10 -5 , N = 74 257, and j3 — 4.51. 
The predictions of the weakly nonlinear theory for these cases were discussed above in conjunction 
with figure 2. 

In the first case, only a Dean disturbance is included in the initial conditions. The weakly 
nonlinear theory suggests that an equilibrium state with B — 7.0 X 10 - ' 5 develops corresponding to a 
maximum streamwise velocity perturbation of 5.84 percent of U . We have initiated the simulation 
with the Dean linear eigenfunction having an initial strength, B — 6.95 X 10 - ' 5 . At t — 3.4845 X 10 4 , 
the magnitude of the instantaneous growth rate has decreased by more than a factor of 100 and the 
simulation is stopped with B — 6.7 X 10 - ' 5 . In this case, the difference between the equilibrium state 
predicted by the weakly nonlinear theory and that obtained from the DNS is less than 5 percent. 

In the second case, only a TS wave with initial amplitude A — 3.5 X 10 -3 is included. The weakly 
nonlinear theory predicts a TS equilibrium state with A — 3.8 X 10 -3 , resulting in a maximum 
streamwise disturbance of 1.94 percent of U . The DNS is stopped at a time t — 1.1367 X 10 4 
corresponding to over 480 TS wave periods. At this time, A — 4.0 X 10 -3 . The difference between 
theory and DNS is 5 percent. 

The remaining two simulations with these parameters include both Dean and TS perturbations in 
the initial conditions. In both cases we have initiated the Dean disturbances with B — 6.95 X 10 - ' 5 . 
This value is slightly less than the equilibrium amplitude that we obtained in case 1. 

When the initial strength of the TS disturbance is A — 2.94 X 10 -4 , the weakly nonlinear 
phase-plane diagram in figure 2 shows that the disturbances are in a region of attraction of the 
saturated Dean vortex. Here the TS wave is expected to decay and the Dean vortex should grow 
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to its equilibrium state. In figure 4 we plot the amplitude of both the Dean and TS disturbances, 
normalized by their respective amplitudes at time t — 0. The long-time prediction of the weakly 
nonlinear theory agrees with the results of the DNS. 



Figure 4. Time evolution of TS and Dean disturbances for R = 629F67, A = 2.189 x 10 5 , N = 74 257, and /l = 4.51. 
The initial TS disturbance has A = 2.7 x 10 -4 ; the initial Dean disturbance has B = 7.0 x 10 -5 . 



Figure 5. Time evolution of TS and Dean disturbances for R = 6291.67, A = 2.189 x 10 ®, K = 74 257, and /l = 4.51. 

The initial TS disturbance has A = 2.7 x 10 -3 ; the initial Dean disturbance has B = 7.0 x 10 -5 . 

For a stronger TS disturbance (A — 2.94 X 1 0 — 3 ) , the disturbances are in a region of attraction 
of the saturated TS wave. Figure 5 shows the amplitude histories for this case. Here it is the TS 
wave that grows toward its equilibrium value and the Dean disturbance that experiences rapid decay. 
Again the weakly nonlinear theory and the DNS agree. 

Similar comparisons are made with R — 5000, A = 2.179 X 10 -5 , N = 10 5 , and f3 — 2.0. The 
predictions of the weakly nonlinear theory for this set of parameters are summarized in figure 3. 
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A numerical simulation initiated with A — 2 X 10 -3 and B — 0 confirmed that the amplitude 
A increased with increasing growth rate, at least until A — 4.8 X 10 -3 when the simulation was 
terminated. 

A Dean vortex alone in this flow is predicted to reach an equilibrium state with B — 7.7 X 10 -4 . 
This value of B corresponds to a maximum streamwise velocity perturbation that is 62.4 percent of 
the undisturbed laminar bulk velocity. Using this as an initial amplitude, an equilibrium state with 
B — 4.3 X 10 -4 is obtained in the direct numerical simulations. The difference between the DNS 
result and that of the weakly nonlinear theory is not unexpected since the disturbance amplitudes 
are so large. In addition, in this case the Landau coefficients &i o and &o,l are quite sensitive to the 
amplitude definition. This warns us to temper our expectations from the theory for these parameter 
values. 

The final simulation performed used initial values of A and B as A — 2 X 10 -3 and B — 6. 5 X 10 -4 . 
In this case, the weakly nonlinear theory predicts that the TS wave will decay, leaving only a Dean 
disturbance in the flow. The time evolutions of the TS and Dean disturbances are plotted in figure 6 
through t — 0.5 X 10 3 . Though the simulation was terminated before a steady-state solution was 
reached, the prediction of the weakly nonlinear theory is qualitatively supported. After the initial 
transient, the TS wave decays rapidly. Considering the large amplitude disturbances in the flow, the 
theory and the simulation agree well. 



Figure 6. Time evolution of TS and Dean disturbances for R = 5000, A = 2.179 x 10 5 , K = 10 5 , and /l = 2.0. The 
initial TS disturbance has A = 2.0 x 10 -3 ; the initial Dean disturbance has B = 6.5 x 10 -4 . 

3.3. Comparison With Previous Theory 

In this subsection we compare our Landau coefficients with those of DHZ for two specific 
cases for which normalized Landau coefficients are available in table 1 of reference 10. Using our 
nondimensionalizations, the Landau constants from the current theory and from the theory of DHZ 
are shown in tables 3 and 4. The parameters of table 3 were studied extensively by Singer and 
Zang (ref. 13) and are discussed above. The quantities ago and 6 q 0 from DHZ were taken directly 
from their linear eigenvalue solver, rather than by perturbing away from the neutral stability curve. 
Except for the coefficient b \ q, which represents the effect of the TS wave on the Dean disturbance, 
the coefficients show satisfactory agreement. For the cases considered here, the variation in b \ q 
with various amplitude definitions is not large enough to explain the discrepancy. Singer and Zang 
(ref. 13) suggested that a yet-unknown error in DHZ led to erroneous values of b \ q. Our work here 
supports that suggestion. 
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Table 3. Comparison of Landau Coefficients for R = 6291.67, A = 2.189 x 10 ®, 
N = 74 257, and /? = 4.508 


Theory 

“0,0 

“1,0 

“0,1 

o 

o 

&1,0 

bo,i 

Current 

1.21 x 10“ 4 

-8.25 

-95000 

3.97 x 10“ 5 

-482 

-7640 

DHZ 

1.20 x 10“ 4 

-8.85 

-101000 

3.97 x 10“ 5 

1.23 x 10 5 

-8470 


Table 4. Comparison of Landau Coefficients for R = 4175, A = 2.179 x 10 5 , 
N = 86 800, and /? = 2.76 


Theory 

“0,0 

“1,0 

“0,1 

O 

o 

h,o 

bo,l 

Current 

1.18 x 10“ 5 

3.97 

-46600 

7.66 x 10“ 6 

-233 

-3340 

DHZ 

1.33 x 10“ 5 

3.97 

-46900 

7.66 x 10“ 6 

1.44 

-3390 


4. Concluding Remarks 

We have presented a weakly nonlinear theory to describe interactions between Tollmien- 
Schlichting (TS) waves and Dean vortices in curved channel flow. The approach used to calculate the 
Landau coefficients is an extension of that developed by Herbert. This approach provides not only 
the Landau coefficients but also the shape change of the fundamental disturbances with amplitude. 
We also discussed the effects of amplitude definition on the Landau coefficients. 

The analysis indicates that the resultant flow can be quite sensitive to the details of the initial 
conditions. By the addition of a Dean vortex of relatively small amplitude, one can change the 
long-time evolution of the flow dramatically. In one case we observed a switching of the equilibrium 
solution from a finite-amplitude TS wave with no Dean vortex to a finite-amplitude Dean vortex 
with no TS wave. A more dramatic change occurred in another case where a TS wave alone in the 
flow would have grown ad infinitum, but the presence of a Dean vortex caused the TS wave to decay 
and the Dean vortex then evolved toward an equilibrium state. 

Direct numerical simulations verify the predictions of the theory. When the amplitudes of the 
disturbances are not too large, the direct numerical simulations and the weakly nonlinear theory 
agree quite well. When the amplitudes of the disturbances are large (on the order of 60 percent of 
the centerline velocity), the simulations and the theory still qualitatively agree. 

A comparison of our results with those of Daudpota, Hall, and Zang (DHZ), who also developed 
a similar weakly nonlinear theory, indicates that DHZ have an error in the computation of one of 
their Landau coefficients. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
January 9, 1992 
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5. Appendix A 

The Left-Hand-Side Operator L„ / 


The elements of the operator L n / are given as 


L 1 ’ 1 

L 1 ’ 2 

L 1,3 

L 1,4 

L 2 ’ 1 

L 2 ’ 2 

L 2,3 

L 2,4 

L 3,1 

L 3,2 

L 3,3 

L 3,4 

L 4,1 

L 4,2 

L 4,3 

L 4,4 


L 


n,l 


itt-Ug + 
r 



r 


( Lhl 

L 1,2 

L 1 ’ 3 

^’ 4 \ 

L 2,1 

L 2 ’ 2 

L 2,3 

L 2,4 

L 3,1 

L 3,2 

L 3,3 

L 3,4 

U 4 - 1 

L 4,2 

L 4,3 

L 4,4 ) 

i (W 

2 n 
2 + - 

2 K 2 + 1 

1 d 


?*2 r dr 


2 m K 
r 2 R 


d?_ 
dr 2 


0 

d 

dr 


d_u l+ u l 

dr r 


2 m K 
r 2 R 


m-u e + i 

r R 



?i 2 K 2 + 1 


1 d 

r dr 


d?_ 

dr 2 


0 

nK 

i — 

r 

0 


0 


r R 



1 d 
r dr- 



illd 

1 d 

r dr 


i — 

r 


ilfS 

0 



6. Appendix B 

The Right-Hand- Side Vectors R s 

The elements of the right-hand-side vectors R s are given explicitly as 

( R l\ 

R] 

R s = 


Equation system s — 5 corresponds to the TS harmonic. Its right-hand-side vector is given by 

r \ = -(.•N«; mo «; ,0 - 0 -v >•+ («5- 0,0 -°i 2 /r- 

Rj = -mu'^y/r - 

Rl = -( + «;’°’ 0 ’°d;(d’ 0 ’ 0 ’ 0 ) 

«5 = 0 

Equation system s = 6 corresponds to the Dean harmonic. Its right-hand-side vector is given by 

«6 = 0 

Equation system s — 7 corresponds to the product of the TS and Dean fundamentals. Its right- 
hand-side vector is given by 


d 


0,0,1,0\ -i,o, 0,0 




d 


+ ± :K’ U ’ U ’ U ) 


1,0,0, On ,, ,0,0, 1,0 


y/r ’ ' " r dr 

V 1 - 11 - 1 - 11 ! 

dr 0 ’ dr 0 ’ 

f?| = i/K ';"'"' 1 '"'; 1 ;,b,"- , '- 1 V- 11 -"-" [j' „i.o.o.o+ i- (u i,o,Q, 0) u o,o,i,o 


dr 


= 0 
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Equation system s — 8 corresponds to the product of the TS fundamental and the complex 
conjugate of the Dean fundamental. Its right-hand-side vector is given by 

III = i /'K'/ 1 ' 11 — l ' ll '( l - ll - ll -"i/r • i l - 11 ,, 1 - 11 - 11 - 11 . iv,,"-"- 1 

4 = - ( i K u ° _ 1 ’% l’ 0 ’°’°)/r + - («J’° + «J’ 0 ’°’°«°- 0 - 1 -°)/r 

dr 9 ‘ dr 0 ’ 

4 = j/K'/ 1 ' 11 — l ' ll '( l - ll -"- ll i/r • i 1 1 111111 


* (1 


0, 0,-1, On ,, ,1,0, 0,0 




1 , 0 , 0 , On ,, 0 , 0 , - 1,0 


f?| = 0 


Equation system s = 9 corresponds to the self-correction of the TS fundamental. Its right-hand- 
side vector is given by 

( 2 u - 1 ’ 0 ’ 0 ’ 0 « 2 ’ 0 ’ 0 ’°+ 2u 0 ’ 1 ’ 0 ’V’ 0 ’°’ ( Vr 

i yA u q uq & ti Q Uq jj i 

dr dr dr dr 


( ■ \ 1,0, 0,0 

(a 1,0 - * 0i , o .) V 


1 , 0 , 0,0 2 , 0 , 0,0 0 , 1 , 0,0 1 , 0 , 0,0 


n>2 _ -nn, -i,u,u,u z,u,u,u 
I lg — — Uq 


Ur 


Ur 


0, / 

.)/*’ 


, — 1 , 0 , 0,0 2 , 0 , 0,0 , 2 , 0 , 0, 0 - 1 , 0 , 0,0 , 0 , 1 , 0,0 1 , 0 , 0,0 , 1 , 0 , 0,0 0 , 1 , 0 , On / 

(Uq u r ’ ’ ’ + Uq u r ’ ’ ’ +Ug u r ’ ’ ’ + Uq u r ’ ’ ’ )/r 




d 


dr v dr 

1,0, 0,0 


dr 


dr 


- ( a l ,o - *01,0.) V 


-f (ur 1 ’ 0 ’ 0 ’ 0 ) « 2 ’°’ 0 ’° + -f (« 2 ’ u ’ u ’ u ) u4^+ ■; i '/'•'•"•"i + „ 


2 , 0 , 0 , On — 1 , 0 , 0,0 


0 , 1 , 0 , On 1 , 0 , 0,0 


1 , 0 , 0 , On 0 , 1 , 0,0 


dr 


dr 


dr 


dr 


( ■ \ 1,0, 0,0 
( a 1 ,0 - *01,0.) V ’ ’ 


^9 = 0 


Equation system s — 10 corresponds to the self-correction of the Dean fundamental. Its right- 
hand-side vector is given by 
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Equation system s — 11 corresponds to the correction of the TS fundamental by the Dean vortex. 
Its right-hand-side vector is given by 
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Equation system s — 12 corresponds to the correction 
Its right-hand-side vector is given by 


of the Dean fundamental by the TS 


LI I LI I LI I 


+ ^4 X0 -°) £( u7 w -°+ «r 1 ' 0 ' 0 ' 0 ] - (ii.o - (*,») 


; 0 , 0 , 1,0 


' dr y r ' r dr r 

D 2 _ — 1,0, 1,0 1,0, 0,0 , 0,0, 1,0 
ri-^2 0 




12- ~>-P \ u e u z' ' ' t u e u z' ' ' + u e ' ' ' ) ~ \ u e *v ' ' 

. - 1 , 0 , 1,0 1 , 0 , 0,0 , 0 , 0 , 1,0 0 , 1 , 0,0 , 0 , 1 , 0,0 0 , 0 , 1,0 , 1 , 0 , 0,0 - 1 , 01,0 , 1 , 0 , 1,0 - 1 , 0 , 0 , 0 \ / 
+ Uq U r ’ ’ ’ + Uq U r ’ ’ ’ + Uq U r ’ ’ ’ + Uq U r ’ ’ ’ +Ug U r ’ ’ ’ )/?* 


Ap,,- 1 -0,0,0, 1,0, 1,0 , — 1 ,°,i ,0, 1,0, 0,0 , ,0,0,1,0s 0,1, 0,0 , , 0,1,0,0s 0,0, 1,0 

dl S U 9 > u r + dr \ u e > u r + dr \ u e > u r + dr \ u e > u r 
. d i 1 , 0 , o,o, -i,o,i,o , d 1 , 0 , 1 , o, ■ - 1 , 0 , 0 , ol n -i , , o,o,i,o 

+ ) U r + ~^S U e ) U r - ( b l,0 - ** 1 , 0 .) u 0 

u : ;., = /Ki ,, -1 - 11 - 11 - 11 ,, l - 11 - 1 - 11 . l( l,0,0,0 u -l,0,l,0_ u l,0,l,0 u -l,0,0,0 )/r 

- i ■'( ,,- | - n - 1 - 11 ,, + U 0 , 0 , 1 , 0 U 0 , 1 , 0,0 + yl-O-l-Oy-l, 0 , 0 , 0 ) 

_ d i — 1 U 1,0,1,0 + d (u _ 1,01,0) u l,0,0,0 + d (u o, 0,1, 0) l( 0,l,0,0 + i_ (u 0, 1,0,0) l( 0, 0,1,0 

LI I LI I LI I LI I 

. d io,o,o, - 1 , 01,0 , d , io,i,o, - 1 , 0 , 0,0 ,, -i, , o,o,i,o 

+ ~^S U Z ) u r + ( U Z ) u r - (.* 1,0 - ** 1 , 0 .) u z 


, 0 , 01,0 


*?12 = 0 
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